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Abstract. Sommerfeld paradox roughly says that mathematically Couette 
shear flow is linearly stable for all Reynolds number, but experimentally it is 
unstable to any size perturbation when the Reynolds number is large enough. 
Our study here focuses upon a sequence of 2D oscillatory shears which are 
the Couette linear shear plus small amplitude and high frequency sinusoidal 
shear perturbations. The sequence of oscillatory shears possesses two intrigu- 
ing properties: (1). in the fluid velocity variables, the sequence approaches the 
Couette linear shear, thus it can be viewed as Couette linear shear plus small 
noises; while in the fluid vorticity variable, the sequence does not approaches 
the Couette linear shear; (2). unlike the Couette linear shear, the sequence 
of oscillatory shears has inviscid linear instability; furthermore, with the se- 
quence of oscillatory shears as potentials, the Orr-Sommerfeld operator has 
unstable eigenvalues when the Reynolds number is large enough, this leads to 
transient nonlinear growth which manifests as transient turbulence as observed 
in experiments. 



1. Introduction 

The most influential paradox in fluids is the d'Alembert paradox saying that 
a body moving through water has no drag as calculated by d'Alembert [3] via 
inviscid theory, while experiments show that there is a substantial drag on the 
body. The paradox splitted the field of fluids into two branches: 1. Hydraulics 
— observing phenomena without mathematical explanation, 2. Theoretical Fluid 
Mechanics — mathematically predicting phenomena that could not be observed. A 
revolutionary development of the boundary layer theory by Ludwig Prandtl in 1904 
resolved the paradox by paying attention to the substantial effect of small viscosity 
in the boundary layer [13) . Prandtl's boundary layer theory laid the foundation of 
modern unified fluid mechanics. 

Sommerfeld paradox has the potential of being the next most influential para- 
dox in fluids. The paradox says that Couette flow is linearly stable for all Reynolds 
numbers as first calculated by Sommerfeld [16j , but experiments show that Couette 
flow is unstable when the Reynolds number is large enough. This paradox lies at the 
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heart of understanding turbulence inside the infinite dimensional phase space. Dy- 
namical system studies on the Navier-Stokes how in an infinite dimensional phase 
space is still at its developing stage [10] . A typical dynamical system study often 
starts from fixed points (steady states) and then pursues their invariant manifolds 
to understand the phase space structures. Certain techniques e.g. Melnikov inte- 
grals can then be used to detect the intersection between invariant manifolds and 
the existence of homoclinic or heteroclinic orbits. In some cases, chaos (turbulence) 
can be rigorously proved to exist using e.g. shadowing techniques [llj . Due to the 
large dimensionality of the problem we are having, numerics is often necessary. 

In recent years, there has been a renaissance in numerical dynamical system 
studies on fluids, especially on the Sommerfeld paradox . For plane Couette flow, 
plane Poiscuillc flow, and pipe Poiscuillc flow, unstable steady states with three 
dimensional spatial patterns have been intensively explored numerically [4] [H] [18] 
[17j [6]. These steady states have a universal "streak-roll-critical layer" coherent 
structure often observed in transient turbulence experiments [7]. The main crit- 
icism on such steady states being responsible for transition from the linear shear 
to turbulence, comes from the fact that in the phase space these steady states are 
quite far away from the linear shear. Experiments show that no matter how small 
the initial perturbation to the linear shear is, a transition always occurs when the 
Reynolds number is large enough. 

Here we study a sequence of 2D oscillatory shears. These oscillatory shears 
are built from single Fourier mode modifications to the linear shear. The sequence 
of 2D oscillatory shears approaches the linear shear in the velocity variable but 
not the vorticity variable. All these oscillatory shears have a linear instability 
with a non- vanishing growth rate when approaching the linear shear. We believe 
that such an instability better explains the initiation of transition from the linear 
shear to turbulence. Of course, our oscillatory shears are linearly unstable to 3D 
perturbation modes too. In [12] . we show that 3D shears in a neighborhood of our 
oscillatory shear is linearly unstable too. 

Explorations on two dimensional steady states turn out to be not successful 
[2] [5] . That is, the counterpart of the 3D upper or lower branch steady state has 
not been found in 2D. On the other hand, numerics shows that transitions still 
occur from the linear shear to turbulence in 2D. Near the 2D oscillatory shears, 
inviscid two dimensional steady states (with a cat's eye structure) can be established 
rigorously |12) . Unsteady viscous cat's eye coherent structures revealed in the 
current study are viscous continuations of the inviscid steady states. 

2. The Sequence of Oscillatory Shears 

Two dimensional plane Couette flow is governed by the Navier-Stokes equations 
with specific boundary conditions, 

(2.1) Ui,t + UjU it j = -p ti + euijj , u iti = ; 

defined in a horizontal channel, where Ui (i = 1,2) are the velocity components 
along x and y directions, p is the pressure, and e is the inverse of the Reynolds 
number e = 1/i?; with the boundary condition 

(2.2) ui(t,x,0) = 0, u 1 (t,x,l) = l, u 2 (t,x,0) = u 2 (t,x, 1) = 0; 

and all the variables are periodic in x. By a theorem of Romanov [14] . the Couette 
linear shear [u\ = y, u 2 = 0) is linearly and nonlinearly stable for any Reynolds 
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number (including infinity). The nonlinear stability is in terms of the L 2 
vorticity. Under the 2D Euler dynamics, in terms of L 2 norm of velocity, all the 
2D shears are nonlinearly unstable by a theorem of Shnirelman [15j . 
In this study, we focus on the sequence of oscillatory shears 

c 

(2.3) m = U(y) = y + - sin(4n7ry), u 2 = 0; 

n 

where c is a constant. One can view this as a single mode of the Fourier series of 
all 2D shears satisfying the above boundary conditions (|2.2j) . 

y + c m sin (my). 

m— 1 

As n — > oo, the oscillatory shears approach the linear shear U (y) — > y. On the other 
hand, in the vorticity variable, the oscillatory shears do not approach the linear 
shear U y = 1 + Acir cos(4n7ry) -/-* 1. Thus in the velocity variables, the oscillatory 
shears can be viewed as the linear shear plus small noises. These oscillatory shears 
are linearly unstable under the 2D Euler flow |12) when 

(2-4) £t-<c<T"- 

In this case, the oscillatory shears bifurcate into steady states of 2D Euler flow 
with cat's eye structures |12) . These claims are in consistency with the Romanov's 
theorem and the Shnirelman's theorem mentioned above. First of all, in the vor- 
ticity variable, the oscillatory shear is not in a small neighborhood of the linear 
shear, while Romanov's nonlinear stability only deals with a small neighborhood 
of the linear shear. This leaves the room for our oscillatory shears to be linearly 
and nonlinearly unstable. In the velocity variable, our oscillatory shears can be 
in a small neighborhood of the linear shear, but Shnirelman's nonlinear instability 
result claims that practically all 2D shears (including the linear shear) are non- 
linearly unstable in the L 2 norm of velocity. Such a nonlinear instability is not a 
result of linear (exponential) instability rather some kind of linear in time growth. 
Such a nonlinear instability does not seem to be responsible for the transition to 
turbulence observed in experiments. The L 2 norm of velocity is more relevant to 
weak solutions which have little physical meaning. We believe that the instability 
of our oscillatory shears with the peculiar feature of approaching the linear shear 
in the velocity variable but not in the vorticity variable, offers a better explanation 
to the initiation of transition from the linear shear to turbulence than the upper 
and lower branch argument. The 3D nature of the transition in experiments can be 
addressed too. Of course, our oscillatory shears are linearly unstable to 3D pertur- 
bations. In fact, 3D shears in a neighborhood of our oscillatory shears are linearly 
unstable too |12j . 

Under the 2D Navier-Stokes flow, these oscillatory shears are not steady, rather 
slowly drifting, 

(2.5) UCt, y) = y + -e"^ 4 ™^) 2 * s in(4mry). 

n 

Nevertheless by simply using the oscillatory shears (|2.3p under condition (|2.4|) as 
the potential, the Orr-Sommcrfeld operator (linear Navier-Stokes operator) has 
unstable eigenvalues which converge to those of linear Euler operator as e — > + 
. A more precise statement of the Sommerfeld paradox is as follows: 
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• Mathematically, the Couette shear is linearly and nonlinearly stable for 
all Reynolds number R, in fact, all the eigenvalues of the Orr-Sommerfeld 
operator satisfy the bound A < —C/R for some constant C [14) . 

• Experimentally, for any R > 360 (where R = j- in our setting pQ), there 
exists a threshold amplitude of perturbations, of order 0(R~^) where 
1 < M < X depends on the type of the perturbations [9] , which leads to 
turbulence. 

Our main idea of investigating the oscillatory shears (|2.3[) is as follows: The lin- 
ear instability mentioned above will lead to a transient nonlinear growth near the 
oscillatory shears (and the Couette shear). Such a growth will manifest itself in 
experiments as transient turbulence. Here the amplitude of the perturbation from 
the Couette shear will be measured by the deviation of the oscillatory shears from 
the Couette shear and the perturbation on top of the oscillatory shears. 

3. Numerical Results 

For all the simulations, we choose as initial conditions by adding random per- 
turbations to the oscillatory shears (|2.3[) . The L 2 norm of the random perturbations 
is 0.01. Notice from (|2.4[) that the parameter c is in the range (0.04,0.08). For all 
the simulations, we choose c = 0.07. In order for the random perturbations not 
to overwhelm the oscillatory part of our oscillatory shear (|2.3p , n needs to be less 
than 7. Here we shall simulate n = 1,2,3. From (|2.4| and ()2.5|) . we see that 
the time T for the oscillatory shear to drift outside its unstable regime (|2.4|) can be 
estimated by 

ce -£(4 n7 r) 2 T = 1J_ 

2 4tt' 

That is, 

e(4:niT)^ 4 

This time scale agrees well with the time scale of the first pulse in the L 2 norm 
evolution of the random perturbations. For example, in the setting of Figure [3][a), 
T « 35 which almost coincides with the duration of the first pulse. In Figure [lj we 
show a few time flashes of the velocity field evolution. One can see clearly the devel- 
opment of streak and wave modes. If we chose 3D random perturbations, we would 
have observed the roll mode too. As mentioned before, the streak-roll-wave modes 
are generic from the Fourier series expansion. The point we want to stress here 
is that the linear instability of our oscillatory shears (and their 3D neighborhood) 
provides the initiation of transition from the linear shear to turbulence. In Figure 
[U we show the velocity L 2 norm deviation of the solution from the linear shear. 
One can see the development of pulses. As mentioned above, the duration of the 
first pulse agrees with the time scale (|3.1|) for the oscillatory shear to drift outside 
its unstable regime. Comparing Figure Q] with Figure [21(a), one can see that the 
development of streak- wave structure in Figure [T] corresponds to the development 
of the first pulse in Figure [Ha). We can measure the growth rate of the uphill of the 
first pulse. For example in Figure HJa), let m be the first minimum, M be the first 
maximum, and At be the duration from the first minimum to the first maximum, 
then the exponential growth rate a is defined by 
. . 1 . M 
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Figure 1. The development of coherent structures and transient 
turbulence with initial condition (|2.3p plus a random perturbation 
(of L 2 norm 0.01) for n = 1, c = 0.07, and i? = 10000. Here the 
linear shear has been subtracted. 



This exponential growth rate measures the exponential growth in time, of the initial 
deviation from the linear shear. In the cases we simulated, the exponential growth 
rate can be as large as a = 0.25. The exponential growth rate here can also be 
regarded as a measure on how fast the flow is initially leaving the linear shear to 
transient turbulence. In Figure [31 we show the velocity L 2 norm deviation of the 
same solution but from the moving frame of the slow drifting of the oscillatory 
shear (|2.5[) . We denote by uq the velocity field given by (|2.5[) . Comparing Figure 
[2] and Figure [3l one can see that the time instant of the peak of the first pulse 
agrees quite well for n = 1, 2. For n — 3, the initial random perturbations almost 
overwhelm the oscillatory part of the oscillatory shear. That is why the agreement 
is not that good. Such an agreement shows that the first pulse is generated by the 
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(c) n = 2,R= 10000 




(d) n = 2,R = 20000 



(e) n = 3, R = 10000 



(f ) n = 3, R = 20000 



Figure 2. The growth of the L 2 norm of the deviation from the 
linear shear with initial condition (|2.3[) plus a random perturbation 
(of L 2 norm 0.01). The growth rate a is defined as the log e of 
quotient (of the first maximum and the first minimum) divided by 
the time spent, (a), a = 0.044, (b). a = 0.055, (c). a = 0.24, (d). 
a = 0.25, (e). a = 0.17, (f). a = 0.18. 



linear instability of our oscillatory shear; and the slow drifting (|2.5p has very little 
effect on the development of the first pulse. We can also measure the exponential 
growth rate of the uphill of the first pulse as for Figure [2] using the same definition 
(|3.2p . Here the exponential growth rate approximates the unstable eigenvalue of 
our oscillatory shear. One can see that all the exponential growth rates in Figure 
|3]are very close to each other. When n is the same and the Reynolds number R is 
different, this fact is expected since the unstable eigenvalue of our oscillatory shear 
approaches its inviscid unstable eigenvalue as R — » oo. Figure[3]also reveals that the 
unstable eigenvalue does not depend on n substantially. This fact is very important. 
It says that as n increases, our oscillatory shear approaches the linear shear, at the 
same time maintaining a non-shrinking growth rate. This implies that no matter 
how small the initial perturbation to the linear shear is, a transition always occurs 
when the Reynolds number is large enough. The invariant manifold argument 
of the 3D steady states cannot offer such a good explanation. The exponential 
growth rate here measures how fast the initial random perturbation grows when 
observed on the moving frame of the slow drifting (|2.5|) toward the linear shear. 
Next we explains why we observe pulses in Figures [2] and [3] The uphill of the 
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Figure 3. The growth of the L 2 norm of the deviation from the 
slow drifting uq (|2.5p with initial condition (|2.3[) plus a random 
perturbation (of L 2 norm 0.01). The growth rate a is defined as 
the log e of quotient (of the first maximum and the first minimum) 
divided by the time spent, (a), a = 0.11, (b). a = 0.12, (c). 
a = 0.13, (d). a = 0.14, (e). a = 0.11, (f). a = 0.13. 



pulse is generated by nonlinear growth induced by the unstable eigenvalues. Since 
our Reynolds number is quite large, the Navier-Stokes equation is "near" the Euler 
equation, and there is a stable eigenvalue corresponding to an unstable eigenvalue 
due to the near conservativeness. The stable eigenvalue corresponds to the downhill 
of the pulse. Concerning other pulses in Figures [2] and El e.g. Figure E)(e)(f), since 
the random perturbations in these cases almost overwhelm the oscillatory part of 
our oscillatory shear, other instability sources play significant roles. The first pulse 
is due to the linear instability of our oscillatory shear. When the first pulse is 
finished, our oscillatory shear has drifted outside its unstable regime (|2.5[) . but the 
new transient state can pick further instability and develop even higher pulses. This 
further illustrates our point that the linear instability of our oscillatory shear serves 
as the initiator for transition to turbulence. 

Finally, it is always tempting to try to apply the tools of dynamical systems 
to characterize turbulence. One promising tool is the Melnikov integral. In [10j . 
we built the Melnikov integral for the 2D Kolmogorov flow (i.e. with periodic 
boundary conditions in both spatial directions and an artificial force) . The Melnikov 
integral there was built from the kinetic energy and enstrophy, and evaluated along 
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Figure 4. The modulation of kinetic energy E{t) and enstrophy 
G(t) in time t when n = 1, R = 10000. 



a heteroclinic orbit. For the Couette flow, we have not found any proper heteroclinic 
orbit. Nevertheless, we can still investigate the modulation of the kinetic energy 
and enstrophy. When e = 0, the boundary condition (|2.2[) reduces to just the 
non-penetrating condition 

U2(t, X, 0) = U2(t, x, 1) = 0. 

A direct calculation shows that the kinetic energy 

E = j {u\ + u\)dxdy 

and the enstrophy 

G = J £l 2 dxdy, where — d x ui — d y u\ 

are invariant under the Euler dynamics. When e > 0, they are no longer invariant, 
but their time derivatives E t and Gt should be small for small e: 

Et — 2e J UiUijjdxdy, Gt = 2e J flfljjdxdy. 

In Figure QJ we plot the evolution of E, G, E t and Gt along the orbit in Figured] 
Due to the high-frequncy oscillation nature of the oscillatory part of our oscillatory 
shear (|2.3|) . the modulation of the enstrophy (characterized by Gt) is about 100 
times of that of the knietic ennergy. This shows that during the initiation of 
transition to turbulence, there has been a massive vorticity production. 
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4. Conclusion 

We introduced a new theory concerning the initiation of transition from the 
linear shear to turbulence in contrast to the invariant manifold argument of 3D 
steady states for the Couette flow. Our theory focuses on high spatial frequency 
single Fourier mode modifications to the linear shear, resulting in a sequence of 
oscillatory shears with the following peculiar features: 

(1) The sequence of oscillatory shears approaches the linear shear in the ve- 
locity variable but not the vorticity variable, 

(2) The sequence of oscillatory shears are linearly unstable. As the sequence 
of oscillatory shears approaches the linear shear, their linear growth rate 
does not shrink to zero. 

Our theory claims that such a linear instability offers a better explanation for 
the initiation of transition from the linear shear to turbulence, than the invariant 
manifold argument of 3D steady states (upper or lower branch) [6]. The main 
criticism of the 3D steady state argument is that these steady states locate quite 
far away from the linear shear in the phase space; while experiments show that no 
matter how small the initial perturbation to the linear shear is, a transition always 
occurs when the Reynolds number is large enough. The numerical simulations in 
this study confirm our theory. The 2D nature of our study is not a limitation. In 
fact, 3D shears in a neighborhood of our oscillatory shears are linearly unstable too 
|12) . and our oscillatory shears are linearly unstable to 3D perturbations also. The 
key point for answering the Sommerfeld paradox is whether or not there is a linear 
instability happening arbitrarily close to the linear shear. 
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